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We use dispersion relations to evaluate the 7Z box contribution to parity-violating electron scat- 
tering in the forward limit arising from the axial-vector coupling at the electron vertex. The calcu- 
lation makes full use of the critical constraints from recent JLab data on electroproduction in the 
resonance region as well as high energy data from HERA. At the kinematics of the Q wca k experi- 
ment, this gives a correction of 0.0047tJJ;8°04 

to the Standard Model value 0.0713(8) of the proton 
weak charge. While the magnitude of the correction is highly significant, the uncertainty is within 
the anticipated experimental uncertainty of ±0.003. 
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Amongst the many methods for searching for physics 
beyond the Standard Model, the verification of the pre- 
dicted evolution of the Weinberg angle from the Z-pole 
to very low energies is currently of great interest. In 
particular, the Q W eak experiment at Jefferson Lab [l[ is 
designed to measure the weak charge of the proton using 
parity-violating elastic electron scattering (PVES) from 
the proton to a higher level of precision than previously 
possible. In combination with constraints from atomic 
parity violation Q, Qweak aims to either discover evi- 
dence for new physics beyond the Standard Model that 
leads to parity violation in electron scattering or raise the 
limit on its mass scale to above 2 TeV, complementing 
direct searches at the LHC @, 3 . 

In PVES the parity-violating asymmetry in the t — > 
and low energy limit is related to the weak charge of the 
proton Q v w (jf: 
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where <Jl(r) is the cross section for left- (right-) hand po- 
larized electrons, Gf is the Fermi constant, and a is the 
fine structure constant. The arrow serves to remind that 
this relation is only realized when radiative corrections 
are properly accounted for, in particular, any residual 
dependence on the electron energy E or the momentum 
transfer squared t. Including electroweak radiative cor- 
rections, the proton weak charge is defined at zero energy 
and momentum transfer as |4| 



Qw = (l + Ap + A e )(l-4sin 2 M0) + A^) 
+ l=W + U ZZ + D 7Z (0) , 



(2) 



where sin 2 9w(0) = 0.23867(16) is the weak mixing angle 
at zero momentum, and the corrections Ap, A e and A' e 
are given in Q and references therein. The contributions 
□lyiv and O zz from the WW and ZZ box diagrams 
can be computed perturbatively, while the jZ interfer- 
ence correction □ 7 z(E0 in addition depends on physics at 
low momentum scales [4|-|6[ • The current best theoretical 
estimate from Ref. Q is Q p w = 0.0713(8). 

In Eq. @ we have explicitly introduced a depen- 
dence of D 7 z(i?) on the electron energy E. The energy- 



dependence of the other radiative corrections in Eq. @ 
is not expected to be important at the O(GeV) ener- 
gies relevant for PVES. The Q v w extracted from A FY in 
Eq. (JT|) will then differ from Q v w in Eq. © by an amount 
U lZ (E) — U lZ (Q) : which we refer to in what follows as 
a correction to at the particular kinematics of the 
electron scattering experiment. In general the term 
has contributions from the vector electron-axial vector 
hadron coupling of the Z boson {d^ z ) and from the axial 
vector electron-vector hadron coupling of the Z (D^), 

a lZ {E) = u^ z {E) + u^ z {E). 

Given that the Q W eak experiment has a precision target 
of 4.2% on Q v w if wc arc to draw meaningful conclu- 
sions in relation to the Standard Model it is crucial that 
all the radiative corrections to PVES be under control at 
a level well below this target. The first studies of the box 
corrections 0-Q suggested that they were indeed under- 
stood to the required precision, with the uncertainty on 
the least constrained, D 7 z(0) term being 0.65%. 

In their seminal early work, Marciano & Sirlin Q com- 
puted the □^(O) correction, which is dominant in atomic 
parity-violation experiments at very low electron ener- 
gies. This correction was further divided into a high- 
momentum contribution to the loop integral, computed 
at the quark level, and a low-momentum contribution, 
computed with the nucleon elastic intermediate state. 
The entire uncertainty on the calculation was taken to 
arise from the low-energy component 

In a stimulating new analysis, Gorchtein and Horowitz 
used forward angle dispersion relations to evaluate the 
additional correction, d^ z (E), which is negligible at the 
low electron energies characteristic of atomic parity vi- 
olation. However, at the O(GeV) energies relevant for 
PVES [l[, this correction was found to be large, with 
an uncertainty potentially capable of jeopardizing the 
interpretation of the Q woa k experiment. Recent model- 
dependent analyses of the low-energy jZ contribution, 
involving only nucleon and A intermediate states, find 
smaller but non-negligible effects 

In this work we revisit this new D^ Z (E) radiative cor- 
rection with a detailed evaluation of the inelastic con- 
tributions, taking full advantage of the wealth of data 
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available in the resonance region and from deep-inelastic 
scattering (DIS) at high energy. We find that with these 
new experimental constraints the total D^ Z (E) contribu- 
tion is quite well under control, constituting a correction 
of 0.0047t8;2gJi, or 6.6±£;g %, to Q p w at the Q wcak kine- 
matics. 

The Q WO ak experiment is designed to measure the 
PVES asymmetry from protons at E = 1.165 GeV elec- 
tron energy and very low momentum transfer squared, 
—t = 0.026 GeV 2 [l|. The purely electromagnetic and 
weak contributions are independent of the electron po- 
larization, so that the asymmetry directly measures the 
parity- violating interference between the photon- and Z- 
cxchangc amplitudes. 

For the e(k) + N(p) elastic scattering process in the 
forward angle limit, the Born amplitudes 
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depend only on the electron and proton convection cur- 
rents. The PV Z exchange amplitude Ai z Y incorporates 
the difference between right- and left-polarized electron 
currents, so that A PV = 2M*M Z V /\M 7 \ 2 in Eq. © 
with g e A = —1/2 the weak axial charge of the electron. 

The correction D jZ in Eq. ^ arises from the replace- 
a^pv 



ment M z v 



M P J, so that 
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Corrections from the interference of M z v with the 
two-photon exchange amplitude under the replacement 



M- 



M. 11 vanish in the t 



limit, and there- 



fore do not affect the asymmetry 

Applying Cauchy's integral theorem, the real part of 
nY z (E) can be obtained from its imaginary part using a 
standard dispersion relation, 
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where P is the principal value. The integration over nega- 
tive energies corresponds to the crossed 7Z box diagram, 
with the vector hadron correction even under E' — > —E'. 
Consequently 3FJe □}f z (0) = 0, justifying the neglect of 
this term at atomic scale electron energies. Following 
Ref. Q we compute only the inelastic contributions to 
the dispersion integral; the elastic component has previ- 
ously been computed to be small @, [8Hl0j . 

From the optical theorem, the imaginary part of PV 
*yZ exchange amplitude can t> e written in terms of 



the cross section for all possible final hadronic states 



2 3m7WP| 



d 3 k' 



( Ana\ 



-2Gr 



{2nf2E k , \Q 2 J \ V2 



(7) 



where q — k — k' is the virtual four-momentum trans- 
fer (with Q 2 = -q 2 ), M(M Z ) is the proton (Z-boson) 
mass, and (neglecting the mass of the electron) Lj£ = 

5 Tr [y^'^igv -3175)75^] is the PV leptonic tensor 
arising from the difference between right and left po- 
larized electrons. The symmetric part of the hadronic 
tensor 

~< z AirM 

x(p\[Jt;(x)J z (0) + J^x)J^0)}\p) (8) 

can be written in terms of the interference electroweak 
Fj z structure functions, 
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Making a change of variables 
d 3 k' 1 
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and evaluating k-p — ME in the rest frame of the proton, 
we find for the g^-dcpcndent part: 
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where s = M(M + 2E) is the total cm. energy squared. 
The structure functions F^ z are functions of the ex- 
changed boson virtuality, Q 2 , and of the invariant mass 
W of the exchanged boson and proton (or alternatively 
of the Bjorkcn variable x = Q 2 /{W 2 ~ M 2 + Q 2 )). 
The lower limit of the W integration is given by the 
mass W„ = M + of the pion production threshold, 
and the upper limit of the Q 2 integration is given by 
Q2 iax = 2ME(l~W 2 /s). Our definitions of the structure 
functions coincide with the standard definitions given by 
the PDG [HI, in which at large Q 2 and W the F^ z struc- 
ture function, for example, can be written (at leading 
order) in terms of parton distributions q and q as 
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with weak vector charges g v = 1/2 — (4/3) sin 9w and 
g v = —1/2 + (2/3) sin 2 0^ for u and d quarks, respec- 

fv 



tivcly. In particular, in the limit 2g v — > e q the interfer- 



ence structure functions F^ z 
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We note that the expression (|TTj) is a factor of 2 larger 
than that quoted in Ref. jzj. Because of the importance 
of this difference we have carried out a number of checks 
of our result. Most importantly, we have verified that 
Eq. (|11|) reproduces the known asymptotic limit for a 
point-like hadron Q, as well as the independently cal- 
culated result for an elastic nucleon intermediate state 
[8l-fl(ij. We also point out that the relation between the 
structure functions and the virtual photon total cross sec- 
tions used in Q omits a factor (1 — x) relative to the 
usual definition, which leads to an overestimate of the 
contribution by 30-40% (and even more in the resonance 
region) . 

In the region of low intermediate state hadronic 
masses, W < 2.5 GeV, inclusive scattering is dominated 
by nucleon resonances. While there is an abundance of 
electroproduction data in the resonance region, there are 
no direct measurements of F^ z ■ For transitions to isospin 
7 = 3/2 states, such as the A resonance, CVC and isospin 
symmetry dictate that the weak isovector transition form 
factors are equal to the electromagnetic ones multiplied 
by (1 + Qyy)- For isospin 7 = 1/2 resonances, which con- 
tain contributions from isovector and isoscalar currents, 
using SU(6) quark model wave functions one can verify 
that for the most prominent 7 = 1/2 states the magni- 
tudes of the Z-boson transition couplings are equal to 
the respective photon couplings to within a few percent. 

Following the analyses of Refs. [l2|-[l4j we fit all of the 
available data, including the latest from SLAC [lj| and 
JLab [r| [I?]], using the isobar model for each Q 2 , tak- 
ing into account the contributions from four resonances: 
P 33 (1232), D 13 (1520), F L5 (1680) and F 37 (1950). The 
background contribution is taken to have the functional 
form (1 — x)p /x, which allows for a smooth transition to 
the large- W region. The fit allows the inclusive transition 
form factors for each of the resonances to be constrained 
accurately up to Q 2 w 3 GeV 2 . At larger Q 2 , where 
the resonance transitions are not as well determined, we 
extrapolate the form factors using an exponential form 
(ljjj . This introduces a relatively small uncertainty, as 
the resonance contributions are strongly suppressed at 
large Q 2 . 

The results of our fit for the F2 structure function are 
shown in Fig. Q] as a function of W for several values 
of Q 2 , together with the uncertainties associated with 
the parameters of the fit. For comparison, we also show 
the results obtained by adding the contributions from 
resonances parametrized in Ref. [l9| and the background 
from Ref. [2(| , which were used in the analysis of Ref. Q ■ 
The resonance parametrization fl9j fixes the fit param- 
eters from data at the real photon point, Q 2 =0. To 
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FIG. 

nance region for fixed Q 2 . The data are from JLab (circles) 
QJ] and SLAC (triangles) []]|. The shaded (yellow) band 
between the solid lines represents the uncertainty on our fit, 
while the dashed lines are obtained by summing the reso- 
nance fit from Ref. and the nonresonant contributions 
from Ref. 



obtain transverse and longitudinal cross sections, the 
Q 2 dependence was inferred in Ref. Q using a simple 
ansatz. Comparison with the data in Fig. [1] shows that 
this parametrization does not adequately reproduce the 
experimental results in the resonance region. 

In the DIS region the interference structure functions 
can be expressed in terms of leading twist parton distri- 
bution functions (PDFs). However, the range of integra- 
tion in Eq. (fTTj) extends to low W and Q 2 , beyond the 
region of validity of a PDF description. To proceed we 
follow Ref. and approximate F^ z by their electromag- 
netic analogs Fx 2 at very small x where the light-quark 
PDFs arc approximately flavor independent. Here the 
structure functions are proportional to a sum over prod- 
ucts of weak and electric charges, which for three flavors 
are approximately equal @, El ■ 

For F2 we use the parametrization from Refs. (2lLl22T[. 
which is motivated by Regge theory and valid at both 
high Q 2 and low Q 2 , 



F 2 (x 7 Q 2 ) = A P x- A (l-x) n+i 
+ Anx^il-xT 



Q 2 



Q 2 + Aj 
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where the first term accounts for the Pomeron contribu- 
tion, while the second arises from an effective Reggeon 
exchange. The parameters, and the Q 2 dependence of 
the exponents n and A, are given given in Ref. (2lj . with 
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FIG. 2: As in Fig. Q] but for higher values of W. The data 
are from SLAC (triangles) QJ], NMC (circles) [H and HI 
(squares) [23] . The solid lines show the fit [2ll . |22| | used in this 
analysis, with the (yellow) bands representing the uncertainty. 
The dashed lines are the results of the GVD color dipole model 
[20l|, w hile the dotted lines represent the MRST leading twist 
fit [23]. 



a couple of minor adjustments to better fit the high-Q 2 
HERA data [H] and the MRST leading twist fit HU {viz., 
the normalization of A increased b y 8% and its Q 2 cut-off 
mass decreased by 5% relative to |2l|)- 

At larger x (x > 0.4) the flavor dependence of the 
PDFs renders the interference function ~ 30-40% smaller 
than F2. The electromagnetic structure functions there- 
fore provide an upper limit on F^ Z . However, to obtain 
a more accurate estimate, we use F^ z = (F^ Z / F 2 ) UT F 2 , 
with F 2 given by Eq. (fT4")) and the leading twist (LT) ratio 
constructed from the MRST PDFs [H. As Figs. [2] and [3] 
demonstrate, this procedure yields a very good descrip- 
tion of the W dependence of the available SLAC, NMC 
and HI data 0,111; H3 m the kinematics most relevant 
for uY, z {E), and also gives an excellent description of 
the ZEUS data up to Q 2 « 90 GeV 2 HfJ. The fit (O in 
the DIS region also agrees well with the MRST parame- 
terization [2J] of F2. In contrast, the generalized vector 
dominance (GVD) color dipole model [2(|, used in the 
calculation of Ref. Q , slightly underestimates the data at 
lower W, but exceeds the other fits above W ~ 100 GeV. 

The contribution from the F\ structure function to 
$sm,oY z (E) is obtained from F2 and the ratio R = 
a L /a T = (1 + 4:M 2 x 2 /Q 2 )F 2 /{2xF 1 ) - 1 of the longi- 
tudinal and transverse cross sections. The latter has 
been measured only over a limited range of x and Q 2 : 
however, its contribution is numerically small, espe- 
cially at large Q 2 and W . We use the parameteri- 
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FIG. 3: As in Fig. H but for Q 2 = 5, 10, 25 and 60 GeV 2 . 



zation R = ciQ 2 (cxp(— C2Q 2 ) + c 3 exp(— C4Q 2 )), with 
C1...4 = {0.014,0.07,41,0.8}, which provides a good de- 
scription of available data and has the correct photo- 
production and high-energy limits. This parametriza- 
tion compares favorably with the earlier SLAC fit [28| . 
which included a mild x dependence, but was restricted 
to Q 2 > 0.3 GeV 2 . 

Performing the dispersion integration in Eq. ([6]), the 
result for 5ReDY z (£ l ) is shown in Fig. S] as a function 
of the incident electron energy E. Although the inte- 
gration in principle involves an infinite range of W and 
Q 2 , in practice we find that around 80% of the value of 
$teny z (E) at the energy relevant to Q WG ak comes from 
energies below 4 GeV, where the Q 2 range extends to 
~ 6 GeV 2 , and W to - 3 GeV. This is fortunate as it 
is precisely in this region that a wealth of very accurate 
data exists from JLab [l6|, [TtJ ■ 

In Fig. [J] the nonresonant contribution is small at low 
energies, but rises logarithmically with increasing E. The 
resonance contribution increases steeply to a maximum 
at E ~ 1 GeV, before falling off like 1/E. The resonant 
and nonresonant contributions to 3?e U^ Z {E) are 0.0026 
and 0.0021, respectively, at the energy relevant for the 
Qwcak experiment, E = 1.165 GeV. We should note, how- 
ever, that this separation is somewhat arbitrary, as only 
the physically meaningful, total cross section enters into 
our fit. 

The combined correction to Q^,- at the Q wea k energy is 
then 0.00471^^4, or 6.6±o | % of the Standard Model 
value 0.0713(8) for Q v w , with the error band obtained 
from the uncertainty in the fit parameters using a vari- 
ational method. In comparison, the correction found in 
Ref. was ~ 0.003. The major difference with our value 
arises from the additional factor 2 in Eq. (fTTj) , which has 



been verified by the independent checks discussed earlier, 
as well as our use of more recent electroproduction data 
and the correct relation between structure functions and 
the virtual photon cross section. The correction is im- 
portant for the interpretation of the Q wea kexperiment, 
given its projected uncertainty of ±0.003 [1|. It is also 
critical to the physical interpretation of the experiment 
which is expected to constrain possible sources of par- 
ity violation from beyond the Standard Model at a mass 
scale of > 2 TeV While the uncertainty in the re- 
sult reported here is satisfactory from the point of view 
of Qweak, it can be further reduced by incorporating the 
new inclusive parity- violating data in the resonance re- 
gion which should be taken soon at JLab [29| . 
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arrow at E — 1.165 GeV indicates the energy of the Q wea k 
experiment. 
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